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We derive the two-loop evolution equations for the Cabibbo-Kobayashi- 



(N 

Qh! Maskawa matrix. We show that to leading order in the mass and CKM hier- 



archies the scaling of the mixings |K(fep, [Kbp, |V^dP, l^sP ^^^ of the rephase- 
invariant CP-violating parameter J is universal to all orders in perturbation 
theory. In leading order the other CKM elements do not scale. Imposing the 
constraint Aj, = A^ at the GUT scale determines the CKM scaling factor to be 
~ 0.58 in the MSSM. 



The weak interaction quark eigenstates and the quark mass eigenstates differ in the 
Standard Model as described by the Cabibbo-Kobayashi-Maskawa (CKM) matrix. In this 
paper we show that the scahng of the CKM matrix follows a universal pattern to leading 
order in the mass and CKM hierarchies; namely, the CKM mixing elements that involve the 
third generation and CP violation scale together, while the other components of the CKM 
matrix do not scale to leading order. This makes it much simpler to consider the form of the 
quark mixings at any other scale, in particular at the scale of a grand unified theory (GUT). 
The common scaling is a model-independent feature of the evolution, but the amount of 
scaling can vary between theories. 

The Yukawa matrices U and D can be diagonalized by biunitary transformations 
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The CKM matrix is then given by 



u a 
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The Yukawa matrices evolve with energy scale as determined by renormalization group equa- 
tions (RGE). This in turn determines an evolution equation for the "running" CKM matrix 
V(/i). 

The renormalization group scaling to leading order in the mass and CKM hierarchies can 
be represented schematically in the following way: 
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where the scale /i is the range rrit < fi < Mq with Mq the GUT scale. The CP- violating 
rephase invariant parameter J |l[] also scales as J{Mq) = S{fi)J{fi) to leading order. We 
have defined our scaling factors to be unity at the GUT scale, but one could equally well 
choose any convenient scale. 

The two light generation quark and lepton Yukawa couplings evolve in a common man- 
ner determined by the gauge couplings and traces of the Yukawa matrices, while the third 
generation Yukawa couplings receive additional Yukawa contributions. This implies that the 
ratios A„/Ac, Xd/K, Ae/A^ do not evolve. The scaling pattern in Eq. (J^) violates unitarity 
of V, but only at subleading order. For example the relation \Vud\'^ + iKisP + iKifeP = 1 is 
violated by terms that are neglected to leading order in the evolution of iKidP and |\4sP- 
The elements |Kidp and |\4sP must evolve to subleading order to preserve unitarity. A 
practical strategy is to evolve the small mixings X = |\4bp, Y = |\4sP, Z = iKbp, and J 
which completely determine the other entries in the CKM matrix. 

In terms of t = IhIh/Mq) the two-loop RGEs can be written as 
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where the coefficients Xi, a,, etc. depend upon the particle content of the theory and are func- 
tions of the gauge and Yukawa couphngs, i.e. a, = ai{gf, (7I, g^, Tr[UU^], TrfDD''"], Tr[EE''']) 
and Higgs quartic couphngs. The coefficients Xj do not enter into the running of the CKM 
matrix but do influence the diagonal quark Yukawa evolution; only terms involving a factor 
of DD^ can rotate the U matrix, and only terms with a factor of UU^ can rotate the D 
matrix. In the minimal supersymmetric standard model (MSSM) the other coefficients are 
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In the Standard Model they are given by 
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where 



Y2{S) = Tr[3UU^' + SDD^' + EE+] . 
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The coefficients Xi can be found in Refs. [|,|I. 

Following Ma, Pakvasa, Sasaki and Babu [^^ we find the CKM evolution equation 
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Here «,j, /c = u,c,t, . . .; a, P,^ = d,s,b, . . . We henceforth restrict our considerations to 
the three-generation case. Defining the four independent quantities X = |V^6p, Y = |V^usp, 
Z = iKfep, and the parameter J = lmVy^j^Vcs^*gV*d which can completely specify a unitary 
CKM matrix, the other elements are given by P] 
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where 



K = \XYZ{\ -X-Y){\-X-Z)- J2(l - X) 
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The full evolution equations for X, Y, Z and J are given in the appendix. Keeping only 
the leading terms in the mass (Ac/A(, A„/Ac, A^/A;,, A^/A^ << 1) and CKM {X,Z,J « 1) 
hierarchies, these equations simplify considerably and a universal scaling is found 
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where W^ = \V,,\\ \V^,\\ \Vts\\ \Vu\\ J and 
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where W2 = |VusP, iKdP, iVtbl"^, \Vcs\^, iKidP- One does not need the mixing between the 
first two generations to be small {Y << 1) which makes the universality an especially 
good approximation. To leading order it is only necessary to include the third generation 
Yukawa couplings in A^ and A^. Notice that Eqs. (0)-(^) violate unitarity of V, but only 
at subleading order. The solution of Eq. ( p4D is 
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where S* is a scaling factor defined by 
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This reduces [^^ to the scaling factor y'^{fi) in the one-loop semianalytic treatment (ne- 
glecting Afo and A,-), with 
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The general behavior of S{fi) is determined by the sign of a^ (and perhaps also a„ in models 
where tan/5 is large). In the Standard Model the scaling factors are greater than one since 
the one-loop coefficients a„ and a^ are negative. 

One might naively have expected there to be contributions to the scaling of the Wi that 
are not proportional to Wi ; for example a contribution to the running of | Vub \ ^ of the form 
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AclKfel o^ the right hand side of Eq. (|3^) can be of the same order as the contribution 
'^t l^feP- W^e conclude that no such terms arise. 

We find the following RGEs for the Yukawa couplings 
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where a = e, fi, t. Including only the third generation in the sums, these equations reduce 
to the leading order expressions for At, A^, A,- yielding 
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respectively. For the first and second generations the corresponding expressions are 
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In Figure 1 we show contours of constant S{mt) in the MSSM versus the values of the 
Yukawa couphngs A^ and A;, at scale rrit and also at the GUT scale. The contribution to the 
scaling from At can traded off against the contribution for Ab as indicated by Eq. (0). These 
contours are shown versus rrit and tan/5 in Figure 2. The mb{mb) = 4:A GeV contour, which 
requires Ab(Mg,) = \t-{Mq), is plotted as well. The evolution equation for R{,/r = Aft/A,- at 
one-loop in the MSSM is given by 



E d.9^ + A? + 3X1 - 3A^) . (48) 



dRb/r _ Rb/r 
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where di = (—4/3,0, 16/3). For small tan/3 the bottom and tau Yukawa couplings can be 
neglected, and the scaling of Rb/r factorizes into scaling due to the gauge couplings and the 
scaling factor S due to the top Yukawa. Given a fixed gauge sector scaling, the nib contours 
and the contours of constant S coincide for small tan/3. Note that rUb — 4.4 GeV implies 
S{mt) ^ 0.58. 

The numerical calculations performed here are similar to those described in Ref. 0. 
The input values are ai(Mz)-^ = 58.89, a2(Mz)"^ = 29.75, and a^iMz) = 0.116 for 
the gauge couplings and mbirnb) = -iA GeV, mcirric) = 1.2 GeV, ms(lGeV) = 0.15 GeV, 
?7irf(lGeV) = 0.008 GeV, mu{lGeV) = 0.005 GeV for the running fermion masses. We take 
the lepton masses to be rrir = 1.7841 GeV, ?7i^ = 0.10566 GeV and rUe = 5.1100 x 10~^ 
GeV. The GUT scale Mq is determined as the scale at which unification of ai and 0^2 is 
achieved. Given an input value for tan (3 the input masses and the gauge couplings determine 
the Yukawa couplings at the scale rrit. We integrate the two-loop RGBs for the gauge and 
Yukawa couplings as well as the evolution equations for X, Y, Z, and J given in the appendix. 
Our results are not sensitive to the values of the first and second generation fermion masses 
or to the input CKM magnitudes |Kfe(mt)| = 0.043, |K6(mt)| = 0.0045, |Ks(mt)| = 0.221, 
J{mt) = 1.95 X 10~^. For experimentally acceptable values of the quark masses and CKM 
matrix elements, the exact scaling as given by the equations in the appendix differ from the 
universal behavior described by Eq. ( |5^ by r^O.1%. 
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A good approximation for evolving the CKM matrix is to use Eq. ( |34D to evolve \Vi 



|V^bp and J and leave |V^s|^ constant as in Eq. (|35D. Then calculate the remaining magnitudes 
|\/,JusingEqs. (0)-(|2|). 

One can show that the universal scaling behavior described by Eqs. (p^)-(p5|) is main- 
tained to all orders in perturbation theory. However the quantitative effects of the three-loop 
contribution are generally smaller than the sub-leading contributions in the mass and CKM 
hierarchies. 

We now give a sketch of a proof of the universal behavior at an arbitrary order in per- 
turbation theory. A higher order contribution will be of the general form 

1 
(167r2)'? 

where q>m + n + ... + o + p represents the loop order. There is an analogous contribution 
to dD/dt. The exponents m and p could be zero. The coefficient f"^^---°p is calculable 
in perturbation theory but can be obtained only with tedious effort; it is a function of 
the gauge couplings gt and the sum of the eigenvalues of the Yukawa couplings matrices, 
Tr[UU^], Tr[DD"''], Tr[EE"''], and possibly other couplings like the quartic Higgs coupling 
in the standard model. The term in Eq. (H^) generates a new contribution to Eq. (E4|), 
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Then to leading order, \Vtb\ 
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which has the same form as Eq. (0). Consequently 
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A similar argument applies to the cases m = and/or p = 0. 

In summary we have shown that there is a universal scaling pattern in the evolution of the 
CKM matrix when only the leading order terms in the mass and CKM hierarchies are kept. 
This is a very good approximation given the observed hierarchy of the quark masses and 
CKM matrix elements. This scaling behavior persists to all orders in perturbation theory. 
Imposing a GUT scale constraint \[,{Mq) = At-(Mq) constrains the amount of scaling. For 
mbirrih) = 4.4 GeV the scaling factor of the CKM matrix is S* ~ 0.58. 

After this work was completed we learned that similar conclusions about the scaling of 
the CKM elements have been obtained at the one- loop level by Babu and Shafi [0. 
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APPENDIX 

The evolution equations for jViap can be derived from Eq. ( p4D using the substitutions 
in Eqs. (|27|)-(p2D, as performed by Babu P at the one-loop level. At the two-loop level one 
obtains 
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Replacing Aj ^ Aj, A^ ^ A^ and omitting the second contributions proportional to (IGvr 
one recovers the one-loop results of Babu [^ (our definitions of X and Z differ from Ref. ^). 
The two-loop equations have the same overall structure as the one-loop equations because 
both contain the same factor VipYji^Vja in Eq. (13^) . Eqs. (|53D-(^6D, together with the evolu- 
tion equations for the gauge couplings Qi and the Yukawa couplings Aj, form a coupled set 
of differential equations that can be solved numerically. In their full form these equations 



together with Eqs. (p^)-(|3^) preserve the unitarity of the CKM matrix to all orders in the 
hierarchy. 
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FIGURES 



Fig. 1. Contours of constant S^^'^{mt) in the MSSM are shown for values of \t and A^ at (a) 
11 = rrit and (b) /i = Mq. We have taken a^{Mz) = 0.116. 

Fig. 2. Contours of constant S^^'^{mt) in the MSSM are shown in the rrit, tanP plane for 
a^lMz) = 0.116. The dashed line is the mbimb) = 4.4 GeV contour obtained from the GUT 
scale condition Xi,{Mq) = XriM^). The X marks the spot at which XiIMq) = \i){Mq) = 
XtIMq) for this rrih contour. In the small tan/3 region a linear relationship exists between 
rrit and sin/? for each contour. 
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